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It was pointed out recently that supersymmetry can gener- 
ate flavor-changing gluonic dipole operators with sufficiently 
large coefficients to dominate the observed value of e'/e. We 
point out that the same operators contribute to direct CP 
violation in hyperon decay and can generate a CP violating 
asymmetry A(A'L) in the range probed by the current E871 
experiment. Interes tingly, m odels that naturally reproduce 
the relation A = */ md/m s do not generate e'/e but could 
lead to an A(A°_) of O(10" 3 ). 

PACS numbers: 14.20 J, 12.60 J and 11.30 E 



The origin of CP violation remains one of the out- 
standing problems in particle physics. Until recently the 
only observation of CP violation was in the neutral kaon 
mixing, with a value of e » 2.27 x 1CT 3 exp(i7r/4) mM. 
The KTeV and NA48 collaborations have now reported 
observations of direct CP violation in the neutral kaon 
decay amplitudes || , with the world average value being 
Re(e'/e) = (21.2 ±4.6) x 1CT 4 f§. 

Although this result is not inconsistent with the 
standard model prediction, it can be used to constrain 
other models of CP violation |)|-|7j]. In particular, it 
has been found that there can be large supersymmetric 
contributions to e'/e Depending on which new 

contributions are large, there are different consequences 
for other processes such as rare kaon decays || and 
hyperon decays. 

In this letter we concentrate on the supersymmetric 
scenario in which the gluonic dipole operators can have 
large coefficients. In this case there are potentially large 
contributions to both e'/e [|| and to the CP violating 
asymmetry A(A°_) in hyperon non-leptonic decays. 

Experiment E871 at Fermilab is expected to reach a 
sensitivity of 2x 1CT 4 for the observable (A(A°_) + A(EZ)) 
H . The CP violating asymmetry A(A°_) compares the 
decay parameter a from the reaction A — ► p?r~ to 
the corresponding parameter a in — > pn whereas 
A(H~) is the corresponding asymmetry for the mode 
S~ — > A°7r~. These asymmetries have a very simple 
form when one neglects the small AI = 3/2 amplitude, 
for example |l(J, 



where Si — 6°, 8n = —1.1° are the final state irN 
interaction phases for S and P wave amplitudes with 
/ = 1/2, respectively |nj]. (f> SiP are the corresponding 
CP violating weak phases. Recent calculations suggest 
that the strong scattering phases in the A°ir final state 
of the S decay are small fl^] , and, therefore, the current 
theoretical prejudice is that |yl(A° )| will dominate the 
measurement. The standard model prediction for this 
quantity is around 3 x 10 -5 , albeit with large uncertainty 
]lO|,[L3]]. This suggests that a non-zero measurement by 
E871 will be an indication for new physics. 

A model independent study of new CP violating 
interactions has shown that A(A°_) could be ten times 
larger than in the standard model and within reach of 
E871 |l4j . A particular example of an operator in which 
A(A°_) can be this large is precisely the gluonic dipole 
operator fH] . The results of E871, therefore, can have a 
direct impact on supersymmetric models. 

The short distance effective Hamiltonian for the glu- 
onic dipole operator of interest is, 



H eff - C, 
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m s da^G^t a (l+j 5 )s 
; m s da uv G^t a (l- l5 ) s + h.c, (2) 



where Tr(t a t b ) = S ab /2, and the Wilson coefficients C g 
and C g that occur in supersymmetry can be found in the 
literature [15[, they are 



C a 



m g m s 



m g m s 



(3) 



A(A°_ 



tan(<5ii - Si)sm.(<f>p - </> s ), (1) 



The parameters 8f 2 characterize the mixing in the mass 
insertion approximation |]l5| , and x = m^/m 2 -, with 
m g , m g being the gluino and average squark masses, 
respectively. The loop function is given by, 

„, s 22 - 20x- 2x 2 + (16x-a; 2 + 9)loga; 
Go(x)=x • (4) 

Ref. H has noted that, in this form, Go(l) = —5/18 and 
the function does not depend strongly on x. The effect 
of QCD corrections is to multiply the Wilson coefficients 

by pi 
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(5) 



To calculate the weak phases we adopt the usual 
procedure of taking the real part of the amplitudes from 
experiment and of using a model for the hadronic matrix 
elements to obtain the imaginary part. We use the MIT 
bag model matrix elements of Ref. [^0 17 to find for the 
weak phases 



= -2.9 x 10 7 GeV 

^—Go(x)Im[(S d 12)Li; 
ilix nig 

= -3.4 x 10 7 GeV 



RL 



B x 



G Q {x)\m[ (Sf 2 ) LR + (6f 2 ) RL )B. 
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(6) 



(7) 



We have introduced the parameters B s and B p to 
quantify the uncertainty in these matrix elements. We 
then find, 



A(A°_) SUSY 



a s {mg) 



a s (500 GeV) 
{{2.QB P ~ l.7B s )Im(6 d 2 ) LR + {2.0B P 



/500 GeV^ G Q (x) 
V rn~ a 



Go(l) 
1.7B s )Im(<5? 2 ) flL ) 



(8) 



The matrix element of the gluonic dipole operator of 
Eq. (Q) between two baryon states is calculated with 
the MIT bag model in Ref. |l7j |, and we assume that 
this result is accurate to within factors of two. The S- 
wave hyperon decay amplitude is then obtained by using 
a soft pion theorem which can have 20 — 30% corrections. 
The P-wave hyperon decay amplitude is obtained by 
considering baryon and kaon pole diagrams. A leading 
order calculation of the dominant, CP conserving, P- 
wave amplitudes in terms of (octet) baryon poles alone 
works reasonably well for A decays. However, additional 
contributions are needed to explain the P-waves in other 
hyperon decays | fl8| , and the first non-leading corrections 
to the A decay amplitude are large. An example of 
an additional contribution is the kaon pole, which in 
Eq. (0) accounts for about 20% of the P-wave phase. 
To reflect these uncertainties in our numerical analysis 
we use 0.5 < B s < 2.0, while allowing B p to vary in the 
range 0.7B S < B p < l.3B s . 

In a general supersymmetric model there are also 
contributions to the imaginary parts of the Wilson coef- 
ficients of four-quark operators. Of these, the dominant 
contribution to the CP asymmetry in hyperon decays 
(within the standard model) is due to Oq |13|. We 



have checked numerically, that SUSY contributions to 
C§ (as well as to 6*3,4,5,7) are much smaller than those in 
Eq. (||), for a parameter range similar to that considered 
in Ref. ||. 



Although the asymmetry A(A°_) is due to the same 
\AS\ = 1 interaction responsible for e'/e, the two observ- 
ables are qualitatively different. For e'/e, both the AI = 
1/2 and the AI = 3/2 amplitudes are equally important, 
whereas for A(A°_) only the AI = 1/2 amplitude is 
important. In this case, the interference necessary for 
CP violation takes place between S and P waves within 
the AI =1/2 transition. This sensitivity to differences 
between S and P waves accounts for the different coef- 
ficients multiplying {Sf 2 )LR and {5f 2 )R.L respectively in 
Eq. (^). For this same reason, supersymmetric scenarios 
in which e' is enhanced through AI = 3/2 operators ||||] 
do not enhance A(A°_ ). 

In order to quantify A(A°_) in supersymmetric models 
where the operators in Eq. (j2j) have large coefficients, we 
compare Eq. (|J) with their contributions to e'/e ||, 



a s (rrig) 
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m s (m c ) + m d (m c 



)g^) Bg 

58lm[(S d 12 ) L R-(S d 12 ) RL ) . (9) 



To obtain this expression, Ref. ||] uses the K — > rnr 
matrix element from a chiral quark model calculation 
in Ref. jljj and uses the parameter Bq to quantify the 
hadronic uncertainty. We use the range 0.5 < Bq < 2 
motivated by the bag model result of Ref. |2(J and 
the dimensional analysis estimate of Ref. H]. It is 
interesting to note that (e' /c)susy depends on the same 
combination of the mass insertion parameters as the 
weak phase 4> s in Eq. (||). We require (e'/e)susY to 
be equal to the observed value (i.e., e'/e dominated 
by supersymmetry) or less (i.e., e'/e dominated by the 
standard model). 

Comparing Eqs. (||) and (||) one sees that e'/e and 
A(A°_) are proportional to different combinations of the 
coefficients (S d 2 )LR and {5f 2 )RL- For this reason one 
cannot determine the allowed range for A(A°_) solely in 
terms of e'/e. In what follows, we consider the three 
cases: a) Im(<5f 2 )flL = 0, b) lm(5f 2 )LR — 0, and c) 
lm(Sf 2 )RL = lm(Sf 2 )LR motivated below. 

It is useful to recall the origin of the mass insertion 
parameters (Sf 2 )RL and {8 d 2 )LR- They are the mismatch 
between the quark mass matrix and left-right mass- 
squared matrix for down-type squarks (we restrict our 
discussion to the first and second generations). In many 
theories of flavor with approximate flavor symmetries, 
the Cabibbo angle originates in the down-quark sector, 
and we find the mass matrix to be of the form 



am s \ 2 m s X 



(10) 



where a and b are 0(1) coefficients and A is the sine of 
the Cabibbo angle. The (2,2) element is nothing but the 
strange quark mass itself (ignoring 0(A 2 ) corrections), 
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and the (1,2) element is fixed by the requirement that 
the Cabibbo angle is reproduced. The down quark mass 
is given by rrid = (a — b)X 2 . A case of a — and b = — 1 
naturally reproduces the phenomenologically successful 
relation A ~ ^Jmd/m s and deserves a special attention. 
This arises if the off-diagonal elements originate in an 
anti-symmetric matrix such as in U(2) model | p2| . The 
A dependence of each element is a consequence of the 
approximate flavor symmetry, but the constants a and b 
cannot be determined by symmetry considerations alone 
and hence are model-dependent. The same approximate 
flavor symmetry constrains the form of the left-right 
mass-squared matrix. Therefore, the left-right mass- 
squared matrix for down-type squarks is 



m Ir = msuSY 



am 
bm s X 



A 2 cm, A 



dm s 



(11) 



where msusY is the typical supersymmetry breaking 
scale which we take to be the same as the down-type 
squark mass, and a, b, c, and d are O(l) numbers and 
can be complex. The U(2) model gives b = — c. 



After diagonalizing the quark mass matrix Eq. (10), 
the left-right mass-squared matrix becomes 

/ {d-bd-b + bd)X 2 (c-d)A \ rio , 

m ^Ym s { {l _ b j )x J )■ (12) 



Unless special relations hold between 0(1) coefficients, 
there remain off-diagonal elements which contribute to 
flavor-changing neutral currents. The mass insertion 
parameters for s — > d transition are defined as 



(<&) 



LR 



m s (c - d) 



msuSY 



(6 d 12 ) RL = «; 



LR 



m s (b — bd)* 
msusY 



(13) 



It is amusing that the size of the mass insertion param- 
eters given here generates e' according to Eq. (||) at the 
observed level for msusY ~ 500 GeV and a phase of 
0(1). 

The case a), of lm(Sf 2 )LR ^ and Im(Sf 2 ) R L = 0, 
corresponds to the choice a = 1, b = in the quark mass 
matrix Eq. (pT0[ ) and its counter part in the squark mass 
matrix Eq. ( jll| ) b = is also likely to be zero in this case. 
We still expect c, d to be 0(1) and this case is the most 
conservative one. The case b) is the other possible limit 
where c — d happens to have a negligible imaginary part. 
Im(6 — bd) can still generate an interesting contribution 
to e', while A(A°_) can be much larger in this case. 
Finally, the case c) Im(Sf 2 ) R L — Ini(<5f 2 )i« is motivated 
by the phenomenological relation A ~ y md/m s and 
hence a = 0, b = — 1. The anti-symmetry in Md could 
imply the anti-symmetry in M^ d R: and hence b = —c. 
This is indeed what happens in the U(2) model of flavor 
p2j. In this case, lm(Sf 2 ) LR = m s Im(-6 - d)/m SUS Y, 



while lm(5 d 2 )nL = m s Im(6 + d)* /msusY = m s Im(— b — 
d)/msusY = Im(<5f 2 )i-R- Therefore, there is no parity 
violation in the CP-violating part of the operators and 
hence the contribution to e' identically vanishes p3| . In 
this case, the only constraint on the size of A(A°_) comes 
from e as we will discuss below. 

The operators in Eq. ^ also contribute to e through 
long distance effects and we must check that this con- 
tribution is not too large. The simplest long distance 
contributions arise from 7r°, rj and r) 1 poles as noted in 
Ref. p3. They yield, 



SUSY 



mi 



y/2m K Am mj 
Im((ir \H e ff\K ))(n \H SM \K ) K . 



(14) 



In this expression Am is the — Kg mass difference and 
(tt°\Hsm\K°) w 2.6 x 10~ 8 GeV 2 is extracted from K -> 
7T7T data. We get the matrix clement (7r° |7Y ej ff \K°) using 
the MIT bag model result p7|] . Finally, k quantifies the 
contributions of the different poles, k = 1 corresponding 
to the pion pole. In the model of Ref. [^5| k ~ 0.2 whereas 
the contribution of the 77' alone gives k ~ —0.9 JlOj] . Wc 
use 0.2 < \k\ < 1.0 and demand that this long distance 
contribution to e, 



'SUSY 



a s (500 GeV) J \ 



31 ( 500 GeV A k G (x) 
02 G (l) 



6.4 Im^ 2 ) Lfl 



12) RL 



(15) 



be smaller than 2.3 x 10 3 . This leads to the constraint 
\A(A°_ )| < 7.3 x lQ-^Bp. Note that we allowed the range 
0.35 < B p < 2.6, and hence \A(A°_ )| can be O(10~ 3 ); we 
cannot exclude it up to 1.9 x 10 -3 . The constraint on the 
mass insertion parameters from the short-distance effect 
(e.g., box diagrams) is weaker: (l^{5f 2 ) 2 LR ) 1 ^ 2 < 3.7 x 
10- 3 for m~ g = mq = 500 GeV and (Sf 2 ) LR = (Sf 2 ) RL 

The regions allowed by the three cases discussed above 
are shown in Fig. ^. The case a) with LR contribution 
only is the horizontally-hatched region with the central 
value shown as a solid line, and the case b) with RL 
contribution only is the diagonally-hatched region with 
the central value shown as a solid line. The shaded 
region at the top is excluded by the e constraint, and 
is particularly important for case c) in which there is 
no contribution to e'. It is interesting that the best 
motivated case c) allows a large asymmetry in hyperon 
decay. The vertical band shows the world average for e'/e 
and the region to the right of the band is, therefore, not 
allowed. 

In summary, we have studied the supersymmetric 
contribution to CP violation in hyperon decays from 
gluonic dipole operators. We parameterize the hadronic 
uncertainties with the quantities Be, B 8 , B p and k which 
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we allow to vary in reasonable ranges. We constrain the 
size of the coefficients of the gluonic dipole operators 
with the observed value of e' and predict a range for 
A(A°_) depending on whether the LR or the RL operator 
dominates. We find that the size of A(A°_) can be within 
reach of the E871 experiment. Particularly interesting 
is the scenario c), which explains naturally the relation 
A = y md/m s . This scenario does not generate e', but it 
can lead to an A(A°^) as large as 10~ 3 . 
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FIG. 1. The allowed regions on 

(\{e /e)suSY |, \A(A°_)susy\) parameter space for three cases: 
a) only lm(Sf 2 )LR contribution, which is the conservative 
case (hatched horizontally), b) only Im(5f 2 )RL contribution 
(hatched diagonally), and c) Im(<5i 2 )-L-R — Ini(<5r2)ftL case 
which does not contribute to e and can give a large |A(A^)| 
below the shaded region (or vertically hatched region for 
the central values of the matrix elements). The last case is 
motivated by the relation A = yjm,d/m s . The vertical shaded 
band is the world average Q of e'/e. The region to the right 
of the band is therefore not allowed. 
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